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STIRLING PERMUTATIONS, CYCLE STRUCTURES OF PERMUTATIONS 

AND PERFECT MATCHINGS 

SHI-MEI MA AND YEONG-NAN YEH 


Abstract. In this paper we provide a unified combinatorial approach to establish a connection 
between Stirling permutations, cycle structures of permutations and perfect matchings. The 
main tool of our investigations is MY-sequences. In particular, we discover that the Eulerian 
polynomials have a simple combinatorial interpretation in terms of some statistics on MY- 
sequences. 
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1. Introduction 

Stirling permutations were defined by Gessel and Stanley 0- Let j 2 := {j,j} for j > 1. 
A Stirling permutation of order n is a permutation of the multiset {l 2 ,2 2 ,...,n 2 } such that 
every element between the two occurrences of i are greater than i for each i £ [n], where [n] = 
{1, 2,... , n}. Denote by Q n the set of Stirling permutations of order n. For a = o\a 2 • • • cr 2n € 
Q n , an occurrence of an ascent (resp. a plateau ) is an index i such that a t < ai + \ (resp. 
ai = cl_|_ i), Recently, there is a large literature devoted to Stirling permutations and their 
generalizations. The reader is referred to HI El 0 El m for recent progress on the study of 
statistics on Stirling permutations. 

In this paper, we always assume that Stirling permutations are prepended by 0. That is, we 
identify an n-Striling permutation ot<T 2 • • • a 2n with the word (Jo (J i (T 2 ■ ■ ■ d 2 n , where do = 0. Let 
a = or d 2 • • • (J 2n £ Q n - We say that an index i £ [2 n — 1] is an ascent plateau if d 9 ;_i < d* = dj+i 
(see m ). Let ap (d) be the number of the ascent plateaus of a. For example, ap (221133) = 2. 
We define 

N n (x)= J2 ^ P(<T) - (!) 

cr&Qn 

Analyzing the placement of 2 copies of (n + 1), it is easy to deduce that the polynomials N n (x) 
satisfy the recurrence relation 

N n+ i(x ) = (2 n + l)xN n {x) + 2x(l - x)N' n (x) 

with the initial value Nq(x) = 1. The exponential generating function for N n (x) is given as 
follows (see m Section 5]): 

N (X,Z) = \/ 1 
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The first few of N n (x ) are 

Ni(x) = x , N 2 {x ) = 2x + x 2 , Ns(x) = 4x + 10x 2 + x 3 , N^(x) = 8x + 60x 2 + 36x 3 + x 4 . 


Let & n denote the permutation group on the set [n] and it = tt( 1 )tt( 2 ) ■ ■ ■ ir(n) € & n . An 
excedance in it is an index i such that ir(i) > i. Let exc (-/r) denote the number of excedances in 
7 r. The classical Eulerian polynomials A n (x) are defined by 

Aq(x) = 1, A n (x) = x exc{7T) for n > 1, 
it es„ 

and have been extensively investigated (see [H|5]:.7j for instance). In |5], Foata and Schiitzenberger 
introduced a g-analog of the Eulerian polynomials defined by 


A n (x-,q ) = J2 x exc(7r) g cycW . 

Tr£&n 

where eye (7r) is the number of cycles in it. Brenti mm further studied g-Eulerian polynomi¬ 
als and established the link with (/-symmetric functions arising from plethysm. In particular, 
Brenti [3j, Proposition 7.3] obtained the exponential generating function for A n (x\q): 

1 — x 


1 + ^2A n (x-,q) 


n>l 


e z(x-l) _ x 

For k > 1, the 1 /k-Eulerian polynomials An\x) are defined by 

E4‘>wx 


n> 0 


1 — X 


gkz(X— 1) _ rjQ 


(3) 


Let e = (ei,e 2 , ■ • ■, e n ) G Z n . Let I U) k = {e|0 < e* < (i — 1)At}, which known as the set of 
n-dimensional k-inversion sequences (see [15]). The number of ascents of e is defined by 


asc (e) = ff <. i : 1 < i < n — 11 — 


{i — 1 )k + 1 

Recently, Savage and Viswanathan m discovered that 


< 


e i +1 1 

ik + 1 J 


From ([2]) and (]3]), we get 


Hence 


AW(x) = Y ^ Sc(e) = k n A n (x- 1 Ik). 

€■^1n,k 

4 2) (X) = X U N n Q) . 

N n (x) = Y /- excW 2 n - cyc W. 

7rG6n 


Let N n (x) = x n N n (^). It follows from ([3]) and (]4]) that 


(4) 

(5) 


2 n A n (x) = Y[l)ftk(x)N n - k (x). 

k =0 y ' 

A perfect matching of [2n] is a set partition of [2n] with blocks (disjoint nonempty subsets) 
of size exactly 2. Let M. 2 n be the set of matchings of [2n], and let M G Ad 2 n- The standard 
form of M is a list of blocks (H, ji)/(* 2 >J 2 )/• ■ ■ /{imjn) such that i r < j r for all 1 < r < n 
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and 1 = ii < *2 < ■ ■ ■ < *n- Throughout this paper we always write M in standard form. It is 
well known that M can be regarded as a fixed-point-free involution on [2n] and also as a Brauer 
diagram on [2n] (see [TO] for instance). Let so(M) be the number of blocks of M with odd 
smaller entries. It is well known that 


N n (x)= * S ° (M) > 

M£M 2 n 

which has been studied in mum- In particular, from mi Theorem 9], we have 


( 6 ) 


2 n xA n (x) = ^2 (^j N k{x)N n - k (x). 


Combining JT]), © and ©, we get 


E 


x a P (<?) 


X n-exc (-tr) 2 n-cyc (ir) _ 


E 


x so (M) 


(7) 


oeQn 7r€6„ M£M2n 

It is natural to consider the following question: Is there existing a unified combinatorial approach 
to prove (J7|>? The main object of this paper is to provide a solution to this problem. The paper 
is organized as follows. In Section [21 we present the main results, and recall some definitions 
that will be used throughout the rest of this work. In other Sections, we give combinatorial 
proofs of these main results. 


2. Definitions and main results 

Let N n (x ) = Y'/u—i N{n , k)x k . Hence N(n, k ) is the number of perfect matchings in Ad 2 n with 
the restriction that only k matching pairs have odd smaller entries. We partition the blocks of 
M into four subsets: 

OO(M) = {(a, b) G M | a and b are both odd}, 

OE(M) = {(a, b) E M | a is odd and b is even}, 

EO(M) = {(a, b) E M | a is even and b is odd}, 

EE(M) = {(a, b) E M | a and b are both even}. 

Let 00 (M) = #00(M),oe(M) = #OE(M),eo(M) = #EO(M),ee(M) = #EE(M). It is evident 
that so (M) = 00 (M) + oe (M). Here we present a further characterization of so (M). 

Proposition 1. We have so (M) = oe (M) +ee(M), and the numbers N(n,k) satisfy the recur¬ 
rence relation 

N(n + 1, k) = 2kN(n, k) + (2 n — 2 k + 3)N(n, k — 1). (8) 

for n,k > 1, where N( 1,1) = 1 and N( 1, k) = 0 for k > 2 or k < 0. 

Proof. Given a M G Ad 2 n - Note that 

2oo (M) + #oe (M) + #eo (M) = n = 2#ee (M) + #oe (M) + #eo (M). 

Then jfoo (M) = jfee (M). Hence so (M) = oe (M) + ee (M). Therefore, we have 


N(n, k) = #{M e Ad 2n : oe (M) + ee (M) = k}. 


(9) 
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We now prove (HI). Let (a, b ) be a given subset of M. Let (p be the construction of M' £ M 2n + 2 
that replacing (a, b) by (a, 2n + 1)/(b, 2n + 2) or (a, 2n + 2)/(6, 2n +1). We distinguish two cases. 

(ci) If (a, 6) € OE or (a,b) £ EE, then the construction (p does not increasing the number 
odd smallers. Combining Q, this accounts for 2kN(n,k) possibilities. 

( 02 ) If (a, 6) £ 00 or (a, b) € EO, then the construction p> does form a new odd smaller. 
Moreover, we can also append (2n + l,2n + 2) to M. This gives (2n — 2(k — 1) + 
1 )N(n, k — 1) = (2 n — 2k + 3 )N(n, k — 1) possibilities. 

□ 


For k > 1 and i > 0, we define 

P k = {1,2, 3,... , 2k} and N e = {-1, -2, -3,. .. , -21, a}. 

Let Y n = (yi,y 2 ,...,y n ), where yi £ P k U N e for 1 < i < n. In particular, Pi (J = 
{1,2,*},Pl U-^i = {1, 2, — 1, —2,*} and P 2 (JN 0 = {1,2, 3,4,*}. Let pos(Y n ) (resp. neg(T n )) 
be the number of positive (resp. negative) entries of Y n . Denote by star (Y n ) the number of * of 
Y n . Careful consideration of Proposition [1] yields the following definition. 

Definition 2. We call the sequence Y n a MY -sequence of length n if y± = * and y k £ Pi+s k U ^t k 
for 2 < k < n, where s k = neg (Y k _\) + star (Y k _\) — 1 and t k = pos (Y k _\) for k >2. 

Note that s 2 = t 2 = 0. Hence y 2 £ PiljA^o- For example, (*, 1,—1,2,*) is a MY-sequence, 
while (*, 1, —1, —4, 2) is not since y\ < —2. Denote by y n the set of MY-sequences Y n . Note 
that 1 + s k + t k = k — 1 for k > 2. Therefore, we have 

#Yn = (2(1 + s n ) + (1 + 2t n ))#y n -! = (2 n - l)#Yn-i = (2 n - 1)!!. 

We can now present the first main result of this paper. 


Theorem 3. For n > 1, we have 


'y ^ ^.oe (M)+ee (M) 
M £M 2 n 


^.neg (Y„)+star (Y n ) 

Yneyu 


V&Qn 


( 10 ) 


In this paper we always write ir £ & n by its standard cycle decomposition, in which each 
cycle is written with its smallest entry first and the cycles are written in ascending order of their 
smallest entry. 


Definition 4. Let (c\,c 2 ,..., cf) be one cycle of n. We say that Cj is a cycle descent if Cj > c ]+ \ 
for 1 < j < i. 


Let cdes (7r) be the number of cycle descents of tt. For example, for it = ( 1 , 3,4, 2)(5, 7)(6), 
we have cdes (n) = 1. For 7r £ & n . it is clearly that exc (7r) + eye (tt) + cdes (n) = n. Therefore, 
it follows from ([5|) that 

N n (x) = E ^,cyc (7r)+cdes (tt) ^n—cyc (7r) 

Now we present the second main result of this paper. 
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Theorem 5. For n > 1, we have 

x neg(Y„)yStar(Y„) 2; pos(Y„) _ ^ (2x) cdes ^ (2z) eXC ^ . 

Yn&yn. TT&Gn 

Note that neg (Y n ) + star (Y n ) + pos (Y n ) = n. The following corollary is immediate. 
Corollary 6. For n > 1, we have 


( 11 ) 


Mx)= n 


Y„ey n 


^ \ neg (Y n )+pos (Y n ) 


X 


pos (Y n ) 


Equivalently, 


2 n A n (x) = E 2 s tar (Y n ) x pos (Y n ) 
Y n ey n 


In the following sections, we give bijective proofs of the main results. As a consequence, we 
get a desired proof of (J7|). 


3 . Proof of the left equality of (fTIIl) 

In the following discussion, we always put any 2-elements block of M into exactly a parenthesis 
or a square bracket. For (a, b) £ M, if (a, b) £ OE(M) or (a,b) £ EE(M), then we replace (a, b) 
by [a, b\. Otherwise, the parentheses that contain a and b unchanged. For example, we replace 
(1, 3)/(2,4)/(5, 8 )/(6, 7) by (1,3)/[2,4]/[5,8 ]/(6, 7). Set m = 2(1 + s k + t k ). Let M (P 1+Sk U N tk ) 
be the set of perfect matchings of [m\ with exactly 1 + s*. square brackets and t*. parentheses. 

Now we present a bracket-breaking algorithm (BB algorithm for short): 

(51) Note that M{P\ (J ATo) = [1,2]. We take [1,2] as the starting point, which corresponds 
to the first term of a MY-sequence. Since Pi [J Nq = {1,2, a}, we can perform one of the 
following three operations: 

(ci) Using the element 1 of Pi [J Nq to break the square bracket [1, 2] such that [1,2] is 
replaced by (1,3)/ [2,4], 

(02) Using the element 2 of Pi [J Nq to break the square bracket [1, 2] such that [1,2] is 
replaced by [1,4]/(2,3). 

(03) As for *, we append [3,4] right after [1,2]. 

(52) For k > 2, given a M £ M(P±+ Sk [J A/ fc ). We use entries of P\+ Sk U A \ to break some 
square brackets or parentheses of M. Assume that c, d, e and / are positive integers. For 
the entries of Pi +Sfc , we distinguish two cases: 

(ci) Using the element 2 i — 1 to break the i-tli square bracket with the restriction that 
(i) if the 7-tli square bracket with elements 2c — 1 and 2 d, then we replace [2c — 1, 2 d] 
by (2c —l,m + l)/[2 d,m + 2]; (ii) if the +th square bracket with elements 2e and 
2/, then we replace [2e, 2 f] by (2e ,m + l)/[2/, 771 + 2], 

(C2) Using the element 2 i to break the i-th square bracket with the restriction that (i) 
if the z-th square bracket with elements 2c — 1 and 2 d, then we replace [2c — 1,2 d] 
by [2c — 1, 777 , + 2]/(2 d, m + 1); (ii) if the 7-tli bracket with elements 2e and 2 f, then 
we replace [2e, 2/] by [2e, m + 2]/(2/, m + 1). 

For the entries of Nt k , we distinguish three cases: 
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(ci) For 1 < i < tk, we use the element —2 i + 1 to break the i-th parenthesis with the 
restriction that (i) if the i-th parenthesis with elements 2 c — 1 and 2d — 1 , then we 
replace ( 2 c— 1 , 2 d—1) by ( 2 c— 1 , m+l)/[ 2 d—1, m+ 2 ]; (ii) if the i-th parenthesis with 
elements 2 e and 2 / — 1 , then we replace ( 2 e, 2 / — 1 ) by ( 2 e ,m + l)/[ 2 / — 1 ,m + 2 ], 
(C 2 ) For 1 < i < tk, we use the element — 2i to break the i-th parenthesis with the 
restriction that (i) if the i-th parenthesis with elements 2 c — 1 and 2d — 1, then we 
replace ( 2 c — 1, 2d — 1 ) by [ 2 c — 1 , m + 2 ]/ (2d — 1, m + 1 ); (ii) if the i-th parenthesis 
with elements 2 e and 2 / —1, then we replace ( 2 e, 2 / — 1 ) by [ 2 e, m+ 2 ]/[ 2 / — 1 , m+1 ]. 
(C 3 ) As for *, we append [m + 1, m + 2 ] right after M. 

Given a MY-sequence Y n with neg(Y n ) + star (Y n ) = i. Repeat the BB algorithm n-times, 
we can get a unique perfect matching M of [2n] with oe (M) + ee (M) = i. Conversely, given a 
perfect matching M £ M. 2 n with oe (M) + ee (M) = i. If we delete 2n — 1 and 2n, then we can 
find the n-th element of the corresponding MY-sequence. Along the same lines, we can get a 
unique MY-sequence Y n with neg (Y n ) + star (Y n ) = i. Thus the BB algorithm gives a bijective 
proof of the left equality of (HOD . In fact, using the BB algorithm, we give a bijective proof of 
the following result: 

^ ^.oe (M)+ee (M)^eo (M)+oo (M) _ ^.neg (Y„)+star (Y n )ypos (Y n ) 

MeM 2 n Y n ey n 

Example 7. Let Y§ = (★, 1, — 2,4, — 1) and let M = [1, 6]/[2,8]/(3, 9)/ (4, 7)/[5,10]. The corre¬ 
spondence between Y 5 and M is built up as follows: 

l^[l,2]^(l,3)/[2,4]; 

-2 -A (1,3)/[2,4] 4A [1, 6]/(3,5)/[2,4] = [1,6]/[2,4]/(3, 5); 

4 -A [1,6]/[2,4]/(3, 5) [1, 6]/[2,8]/(4, 7)/(3, 5) = [1,6]/[2, 8]/(3,5)/(4, 7); 

-1 -A [1, 6]/[2,8]/(3, 5)/(4, 7) 4A [1,6]/[2, 8]/(3,9)/[5,10]/(4, 7) = [1,6]/[2, 8]/(3,9)/(4, 7)/[5,10]. 

4. Proof of the right equality of (fill 

Set r = 1 + Sk + tk- Denote by Q ( Pi+ Sk U ^t k ) the set of Stirling permutations of order r with 
exactly 1 + Sfc ascent plateaus. In particular, Q (Pi (J No) = {11}, Q (Pi U N±) = {2211,1221} 
and Q (P 2 |J No) = {1122}. Let us first give a definition of labeled Stirling permutations. 

Definition 8. Let a £ Q (P\ +Sk (jN tk ). If i\ < *2 < • • • < h+s k are the ascent plateaus of a, 
then we put the superscript labels 21 — 1 before ii and 21 after it, where 1 < I < 1 + s^. In the 
remaining positions, we put the superscript labels — 1 , —2, ..., —2 tk and * from left to right. 

For example, the labels of 13324421 are given as follows: 

_1 1 1 3 2 3 _2 2 3 4 4 4 _3 2 _4 1*. 

Given a MY-sequence Y n = (yi,y 2 , ■ ■ ■ ,y n )- Now we present a labeled Stirling permutations 
algorithm (LSP algorithm for short): 
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(51) Since Q(P\ (J AT 0 ) = {11}, we take 11 as the start point, which corresponds to y\ = *. 
Since 11 can be labeled as 1 1 2 1 *, we distinguish three cases: if 1 J 2 = 1, then we get 2211 
by inserting 22 to the position with superscript label 1 ; if yi = 2 , then we get 1221 by 
inserting 22 to the position with superscript label 2 ; if y 2 = *, then we we get 1122 by 
inserting 22 to the position with superscript label *. 

(5 2 ) When 2/2 = 1) we take 1 2 2 2 _ 1 1 _ 2 1 * as the start point. Since 1/3 € Pi |J lYi, we distinguish 
five cases: if 2/3 = 1, then we get 332211 by inserting 33 to the position with superscript 
label 1; if 7/3 = 2, then we get 233211 by inserting 33 to the position with superscript 
label 2; if 2/3 = —1, then we get 223311 by inserting 33 to the position with superscript 
label -1: if 1/3 = —2, then we get 221331 by inserting 33 to the position with superscript 
label -2; if 2/3 = *, then we get 221133 by inserting 33 to the position with superscript 
label -k. 

(53) When 2/2 = 2 , we take _ 1 1 1 2 2 2 ” 2 1 * as the start point. Note that 2/3 G P\ (jNi- Along 
the same lines as the five cases of ( S 2 ), we get 133221,123321,331221,122331,122133. 

(5 4 ) When 2/2 = *, we take 1 1 2 1 3 2 4 2* as the start point. Note that 2/3 G P ‘2 U Ad • Similarly, 
we get 331122,133122,113322112332,112233. 

(5 5 ) Repeat the above procedure, we get all labeled Stirling permutations. 

It is straightforward to show that each such labeled Stirling permutation will be obtained 
exactly once. Indeed, given a labeled Stirling permutation of order n with i ascent plateaus, 
we can just read the indices of ascent plateaus. Deleting the pair (2n)(2n) gives the entry y n . 
Along the same lines, we can get a unique MY-sequence Y n with neg(l^) + star (Y n ) = i. In 
conclusion, LSP algorithm gives the desired proof of the right equality of (11011 . 

Example 9. Let Y$ = (*, 2, — 1, 3, —2,2) and let o = 366314455221. The correspondence between 
Yq and a is built up as follows: 

2 -A 1 1 2 1* o 1221; 

-1 -P 2 1 1 2 2 2 _2 1* & 331221; 

3 -A 1 3 2 3 -1 l 3 2 4 2 -2 l* <$■ 33144221; 

-2 -A 1 3 2 3 -1 l 3 4 4 4' 2 2 -3 2 -4 l* <$■ 3314455221; 

2 3 2 3~ 1 1 3 4 4 4 5 5 6 5 _2 2 _3 2 _4 1* 366314455221. 

5. Proof of m 

Let us first give a definition of labeled permutations. 

Definition 10. Let n G (5 n with p excedances. If i\ < *2 < ■ ■ ■ < i p are the excedances, then 
we put superscript labels —k between ik and n(ik), where 1 < k < p. In the remaining positions 
except the first position, we put the superscript labels 1,2,..., n — p from left to right. 
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In the following discussion, we shall attach superscript labels to all permutations in <5 n . For 
example, the labels of (1,3, 7,4, 2)(5)(6,10, 8)(9,11) are given as follows: 

(1 _1 3 ~ 2 7 1 4 2 2 3 )(5 4 )(6 “ 3 10 5 8 6 )(9 “ 4 11 7 ). 

A labeled permutations algorithm (LP algorithm for short) for generating MY-sequences is 
given as follows: 

Take (l 1 ) as the start point, which corresponds to the first term a of a MY-sequence. Let 7 r be 
a labeled permutation in & n . There are n + 1 permutations of 6 n +i can be obtained from 7 r by 
inserting n + 1 to the positions with superscript labels or as a new cycle (n + 1). We distinguish 
three cases: 

(Ci) If we insert the entry n + 1 to the position with superscript label k, then the n-th term 
of the corresponding MY-sequence is y n € {2k — 1, 2k}. 

(C 2 ) If we insert the entry n + 1 to the position with superscript label then the n-th term 
of the corresponding MY-sequence is y n € {— 2£ + 1, —2£}. 

(C 3 ) If we insert the entry n + 1 at the end of 7 r to form a new cycle (n + 1), then the n-th 
term of the corresponding MY-sequence is a. 

In particular, consider 7r = (l 1 ). Note that (1,2) is obtained from (l 1 ) by inserting 2 to the 
position with superscript label 1. Hence 7/2 £ {1)2} is the MY-sequence that corresponds to 
(1,2). Note that (1)(2) is obtained from (l 1 ) by inserting 2 as a new cycle (2). Hence 7/2 = a 
is the MY-sequence that corresponds to (1)(2). Take (l 1 ) as a start point. Repeat the LP 
algorithm n times, it is easy to verify that each MY-sequence of length n will be obtained 
exactly once. Note that 2 cdes ( 7 0+ exc,r is the number of MY-sequences that corresponds to 7r. 
Therefore, the LP algorithm gives a combinatorial proof of (1111) . 


Example 11. Given ir = (1, 3, 5, 2,6)(4). The LP algorithm can be done if you proceed as 
follows: 

(l 1 ) ^ (1,2) ^ 2/2 e {1,2}; 

(1 - 1 2 4 ) —t (1,3,2) 7/3 £ {—1,—2}; 

(1 - 1 3 1 2 2 )^(l,3,2)(4)^7/ 4 e {*}; 

(1 - 1 3 1 2 2 )(4 3 ) -A (1, 3,5, 2)(4) ^ y 5 e {1,2}; 

(1 _1 3 “ 2 5 1 2 2 )(4 3 ) -A (1,3,5,2,6)(4) 4A7/6 € {3,4}. 


Therefore, the corresponding MY -sequences of ir are given as follows: 


(*j 1) — 1)*) 1) 3); (a, 1, — 1, 
(*j 1) ~2, *, 1,3); (a, 1, —2, 
(a, 2, —1, a, 1,3); (a, 2, —1, 


a, 1,4); (a, 1, —1, a, 2,3); (a, 1, —1, a, 2,4); 

a, 1,4); (a, 1, —2, a, 2,3); (a, 1, —2, a, 2,4); 
a, 1, 4); (a, 2, — 1, a, 2,3); (a, 2, —1, a, 2,4); 


(a, 2, —2, a, 1, 3); (a, 2, —2, a, 1,4); (a, 2, —2, a, 2, 3); (a, 2, —2, a, 2,4). 
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